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(Ue Ene SsiS  —CONnSIStS of two parts. First, an 
experimental investigation of a new device called the 
whaistier nozzle was conducted. Experiments were 
conducted in the areas of nozzie efficiency, mass 
entrainnent, and LlLow Visualization. F iow 
visualization showed the presence of a Coanda type jet 
wea lL interaction in the nozzle collar. JileWag ste 
efficiencies indicated that whistling Goud be 
achieved without much greater iosses than the kasic 
axisynmetric jet. Entrainnent tests were inconclusive 
regarding the whistler nozzle performance. Second, 
Supdersonic flow past an oscillating cylindrical shell 
is analyzed using linearized characteristics methods. 
Pressure distributions and generalized aerodynasic 
forces are calculated and presented for various radius 
to length ratios and reduced frequencies. Good 
agreement is oktained in the two-dimensional limiting 
case with previous vork by Platzer, and an early 


Solution of the steady cylindrical case by Zierep. 
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The pheuomenon of entrainnent is a process through which a 
moving fluid increases its mass flow, by "grabbing up" the 


medium it is flowing in. The moving flvid is called the 


meme £C1UIG, and Maignt be the primary wate jet an an 
Uap The mediuia might be any other fluid, say oil, 


ejector p 
in the case of the ejector pump. 

The entrainment process has found practical use for 
quite a long time. Mieery ticat descriptions Cx the 


entrainment process exist; however,none 1S general enough to 


descripe ali cases of practical interest. Even tne 
Merartiveiy Samsie case of a steady jet in still air is 


Lather dixficult oecause of the many configurations which 
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a continues, perhaps with more intensity than 
previcus]y dué to the recent interest in new V/sToL 
Pemttgerations ehich would employ this process. One such is 
Pirewewavy’S XFV= 124, Which 1S expected to derive a large 
Peer Oretnrust througn the augmentation, or entrainment, 


O© itS primary thrust using wing mounted ejectors. 


The success of concepts such as the augmenter wing 
depends cn many factors, an important one being the rate of 
entrainment cf the primary jet, and another being its 
efficiency. This rate of entrainment is the trate of 
increase cf “et masS flow with axial distance away from the 
jet nozzle. The rate of entrainment and nozzie efficiency 


elemetarecteriactscS Which vary with nozzle configuration. 


As mentioned earlier , there is a large variety of 


steady nezzie configurations. Lobe, slot, and hypermixing 





nozzles are some Which nave drawn recent interest for their 
entrainment capabilities. One should not neglect unsteady 
nozzie configurations, however, as such effects as swirling 
the primary flow, and oscillating the primary flow in a bi- 
Seapto shane: (Lltiadic oscillator), aiso show pitomise in 


their ability to entrain. 


One of the newest nozzie configurations to Surface 1S an 
unsteady nozzle. Perhaps the simplest in construction of 
puuimeuiscteady nozzles, it 1S illustrated in figure 1. The 
whistler nozzle, as it 1s called, produces a loud pure tone 
while in operation. Additionally, it creates a relatively 
large rate of entrainment. It is the whistler nozzle which 


imaeene Subject of this thesis. 
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Figure 1 - WHISTLER NOZZLE CROSS SECTION 
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II. RELATED PAST WORK 


The earliest werk done on the whistler nozzle was by Hill 
Gelemobecne {rererence 1). In 1974, they cited the discovery 
of a new phencmenon, characteristic of a new device which 
they named tke whistler nozzle. Their basic confiacuration, 
Somplliustclated in figure 1, consisted of a convergent 
section, folloved by a section of constant area followed by 
a step change in cross sectional area to another section of 
constant area. The Length of the last section could be 
Varied tc produce loud pure tones of varying frequency, 
eeuougy the excitation of a Standing acoustic ave in the 

a: 


meat Secticn, much S in an organ pipe. Hill and Greene 
observed that the whistler produced a greatly increased 
Mixing rate, and deciaed that this was due to accustic 
Stimudaticn of the jet. Increased mixing rate cf an 
axisymmetric jet due to acoustic stimulation 1S a process 
that had been observed earlier by Cio ana 


Champagne (reference 5). 


ieadohesOne tO | the  HasSice axisymmetric configuration, 
many other configurations were found to produce the 
whistling phenomenon. A few are illustrated in figure 2. 
Common to 11 of these was that the step in cross section 
had to extend completely around the jet. Also noted was tine 
Syuearbent presence of a separation reattachment cycle, 


SCeoUrring at the nozzie lip. 


Ne further study seems to have been performed ccncerning 
the whistler nozzle, save a report in 1975 by Hill and 
Jenkins (ceference 2), in which they reestablished tie 
Wreereting eharacterListics Of the nozzle, and its ability to 


Clit rain more air. 
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Figure 2 - VARIOUS 


13 


i £3 . 
SCuaren wits: Ler 
with rounded 
corners 
.* : ~ 
oe 
Piper cal Exe 
a= Contour Whistler 


ties (Lar CONFIGURATIONS 





Hill and Greene's work has been summarized, since their 
Original report, in comprehensive reviews of the evyjector 
State of tne art by Schum (rererence 3), and by Victs 
(eererence 4). These reports call attention to the whistler 
nozzie's promise in rate of entrainment, and the lack of 


data regarding the whistler thrust efficiency. 


Although no further research directly involving the 
whistler is avaiiable, much work has been done in the area 
of free jets, and entrainment, which might prove valuable in 


a study of the whistler nozzle. 


An opvicusly imoortant area of related work is the 
Seieoltaivity cr free jets to acoustic stimulation. Prebaply 
the best known and most definitive work in this area was 
performed by Crow and Chatipagne (reference 5). By mounting 
a speaker inside the plenua chamber they managed to excite a 
free axisymmetric jet, at various frequencies. The effect 


Mame —>OUnG Yas tO produce a velocity fluctuation in tne 
Wa 


nozzle exit plane. Bete S observed that the turbulent 
structure of the jet was sensitive to the acoustic 
excitation. Through the use of high ~ speed schlieren 


phetography it was observed that under acoustic stimulation 
the turbulent structure of the jet developed a iarge~ scale 
vortex structure, as contrasted witn the small scale vortex 
sheet Sutrounding a normal jet. Attendant to the formation 
of the large-scale vortices, an increase in the mixing rate 


was also cbserved. 


Others, such as Bevilagua and Lykoudis (reference 6), 
have observed the turbulent structure in steady, free jets. 
The presence of a large-scale vortex structure has been 
Observed in Such jets at iow Reynolds numbers. It appears 
that the large-scale structure dissappears as Reynolds 


humber increases. However, according to the findings of 


44 





Crow and Champagne, it may be reexcited through accustic 
Stimulation: Tt may be conjectured fron experimental 
evidence that the entrainment process is a combination of a 
Pmeigersacate "Scooping up! of the ambient fluid, and a smali- 
scaie nibblirg away at the ambient fluid. These entrainment 
processes wceuld corresvond to the large and small scale 


vortex structures around a free jet. 


A further area of related research, perhaps not obvious, 
is the well known Coanda effect. The basic whistler 
configuration sets up a flow situation somewhat akin te that 
meee bastable fluidic elenent, or fluidic oscillatcr. The 
dirference kEetween the two cases is that the fluidic 
oscillatcr cscillates between two positions, while the 
oo te> AUSt CcScillate hetween many circumferential 
MOsitLlons. That such an oscillation takes place seens 
medicated oY the nozzle coniiguration and the separation 
reattachment cycle noted by Hill and Greene 


Reference 10 provides good summary of the mechanisa of 


6F 


the Coanda effect, as weil as its characteristics regarding 
reattachment and sensitivity to sound. Though the following 
comments relate specifically to two- dimensional flows, it 
1s thought that the general behavior and mechanism Should 


Carry over to a three- dimensional situation. 


Figure 3 je wet alee = the Coanda mechanisn in 
two-dimensional flow. As shown, a jet 1s issuing frcn an 
orifice and flowing between two adjacent walls. Because the 
jet wili at some ‘time be nearer to one wall than the cther, 
and because the jet is entraining air, a pressure difference 
across the jet will develop that pulls the jet to the closer 
Wall. This attachment forms a region of circulating f£iow. 
The point of jet attachment is determined by this regicn, in 
Which a balance is struck between flow entrained by the jet 


and flow reinjected by the jet. 


is) 





A Attachment 


Pownm¢ 





REGCiTreulatrtion 


Region 


Figure 3 - COANDA EFFFCT 
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Anything that disrupts this balance causes the attachment 


point of the jet to shift. 


Observation has shown that the attachment point in 
Coanda effect 1S sensitive to jet velocity and sound, within 
a cange of Reynolds numbers. In the range of Reynolds 
number from 200 to 3000, under the effect of increasing jet 
Mem ocrtyY, Oz a SOUnd invut of increasing amplitude, the 
attachment pcint has been observed to move upstrean. pais 
upsStreahk movement appears to approach a limiting position, 


beyond which the point will not move. 


reference 10 relates tne effects of jet velocity and 
outa On the attacanent point to the point of transition 
from laminar to turbuient flow of the jet. Both increasing 
jet velocity and a sound input will move the transition 
ocint cioser to the jet exit. This movement will increase 
the length of turbulent flow across the recirculation 
megiOn. RK turbulent jet entrainus more air than a laminar 
jet, and thus the balance of flow ain the recirculation 
meqeon ) 1S disrupted, with a cesuiting upstream move cf the 


attacnnent pcint. 


1 





Ili. SCOPE OF RESEARCH 


The intent of this research was to provide a clearer 
insight into the workings of the whistler nozzle, and to 
develop an improved understanding of the ways in which one 


may effect controi over the characteristics of a free jet. 


x 


Specific areas of fre 


U 


search were the nozzle thrust 
t2 


efficiency, the nozzle entrainment rate, and visualization 


ay 


n 
Smecone Lip interaction process described py Hill and Greene. 


Accordingly, experiments were conducted to obtain relevant 


{2s 


data < Thrusts were measured on a thrust bed, ass 


entrainment measured in a device patterned after that 
developed by Kicou and Spalding (reference 7), and high 


epeed motion pictures taken of tufts about the nozzle lip. 
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The whistler nozzle used in ail tests is illustrated in 
Prgure 4. It consisted of a small plenun, convergent 
Section, Short constant area section and a moveable collar 
over the constant area Section. The inner whistler diameter 
was 74 inch, while the coliar diameter was 1.5 inches. The 
nozzie configuration was essentially that used by Hill and 
Greene ai their investigations. The characreristic 
dimension used for Reynolds number and nondimensionaliization 
of axial distance from the nozzle was the inner nozzle 


Giameter, 7 inch. 
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Figure CS eee SOLA AND PLENUY 
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Memo eoChiP TION Or EXRERTMENTS 
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Thrust efficiencies were obtained by measuring nezzle 
Ciust. \ comparison was made between measured and 


calculated thrust. 


Test equipment consisted of a thrust bed and the nezzle. 
PQjemenr ist bed was CONTigured to measure horizontal thrust, 
fGeeeeraeOn two tracks, mounted to a sotid foundation. ‘The 
nozzle plenum was clamped to the thrust bed which was in 
Meee coUbca tO the Loundation through a lcad ring, weunted 
With a ioad celi, the device that actually measured thrust. 
Air ¥as supplied to the nozzle from a large tank, connected 
to the plenur through a regulator and a flexible tube which 
led verticaliy away from the plenum, so as not to bias the 
thrust measurements in any way. Calibration was perforned 
With the use of a tray, connected to the bed through a cable 


and pulley arrangement, upon which one placed weigats of 


known magnitude. Peripheral gear requixed included a 
voltage source, digital voltage meter, anda Wheatstone 
Perage. Figure 5 illustrates the set up. Plenum pressures 


were measured using a mercury manometer. 
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Calibration of the thrust bed indicated a certain amount 
of hysteresis in the system, as may be seen in que. 6. 
Not serious at thrusts of about 6 pounds, it di@ introduce 
problems at low thrust levels, and precluded reliable thrust 


meaSurement at low Reynolds numbers. 


Thrust measurements were made for Various COd Mae 
positions, over a range of nozzle exit velocities. Recorded 
data included gage pressure in the nozzle plenum, thrust 
readings from the voltmeter, collar extension, and ambient 


pressure and temperature. 


Measured thrusts were plotted against pressure ratio 
across the nozzle. Figure 7 presents an example curve. 
ieee ticrust was caiculated according to the procedures in 
appendix A. Thrust efficiencies were then obtained and 
plotted against nozzie pressure ratio, exit velocity, and 


Reynoids numher (based on the nozzle diameter). 


oe MASS ENLTRAIN®e NT 


Hfass entrainment waS meaSured directiy, uSing a device 


patterned after one developed by Ricon and Spalding 
(reference 7). It was believed that this device should 
yield more accurate results than an approach using 


integrated velocity protiles, especially in the case of an 
unsteady nozzle. MeaSurementS were made at various axial 


positions, for several collar positions. 


The device used will hereafter be referred to as_ the 
entrainment champer. Tt was deSigned with the aid of data 
from a research report by Peschke (reference 8). Appendix B 


provides a detailed account of the design. 
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@ieeeprinciple cf operation is relatively Simple. Tne 
entrained air of a free jet flows in at a right angle to the 
Jet axis. The entrainmwent chamber simply blocks this inflow 
Witm a circumferential membrane, and seals the ends with a 
peate Sat the nozzle and an orifice at the axial station 
where one wWisnes to measure the total mass flovw. Ideally, 
wok tne entrained air is supplied to the exterior of the 
membrane from a regulated source. As a free jet does not 
Support any axial pressure yradient, the proper amount of 
alr to supply through the membrane is indicated by the 
Som@ePedrance Of any pressure gradient across the orifice. 
Sew tay regulate the membrane flow to obtain this foint. 
Miesemweplimary aspects of the device are illustrated in 


meguce 8. 


Critical design aspects of the entrainment chamber, 
especially the orifice sizes, are treated in appendix B. It 
Meeeolrticient here to simply state that one entrainnernt 
chamber was designed to perform measurement of the total jet 


Mass LlLOW at variouvs axial locations. 


One primary difference between this chamber's gecmetry 
and that of cthers (references 7 and 8) should pe noted. As 
Shown in figtre 8, exit orifices were positioned in the 
porous cylinder at various distances from the nozzle. That 
portion cf the porous cylinder extending apove the orifice 
plate was sealed by a fiberglass cyiinder attached to the 
Seviirce, Thus,in this situation that portion of the jet 
beyond the orifice and inside the porous cylinder exit was 
not ina truly free~jet flow situation. This was as opposed 
to the reference 7 and 8 configurations where tne exit 
Pirbereemconresponded to the exit of the porous cylinder. 


Thus, in those experiments the jet past the orifice was ina 


truly free~jet flow situation. The configuration used in 
this report was adopted in order to facilitate the 
measurement of mass flows at many axial locations. It was 
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hoped that this chamber's geometry would not affect the 
chamber readings greatly, although some effect was to be 


expected. 


Peripheral equipment to the basic entrainment chamber 
ancluded flew regulating, metering, and measuring devices. 
Figure 9 fresents a photograph of the chamber and its 
peripheral gear. Separate air sources, both fed by a large 
tank, were sSuppiied for the nozzie and membrane. Nozzle 
flow was metered through an orifice and centrolled with a 
meeikatobs. Eecause of the much greater mass flow required, 
membrane air waS Supplied through a venturi, followed again 


bya £Liow regulator. Pressures across the oritice wer 


O 


PP eveod wich ah anciined water manometer. Peas o ies ac Loss 


the venturi were measured with a Mercury nanoneter. 
Pressure across the entrainment chamber orifices aas 
measured with an aicnohol nicromanonmeter. Nozzie vienui 


pressure waS measured with a water manometer. 
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Nembrane air temperature was measured just after the 


venturi, uSing a mercury thermometer. 


Entrainment meaSurementS were made at various axial 
Mee@ation=> fOr Various collar positions. Recorded data 
included pressure difference across the venturi and orifice, 
pressure difference across the entrainnent chamoer orifice, 
nozzle plenum pressure, ambient pressure, air temperature in 
the membrane Spo Ly; and COmlar YeSsverom Standard 
procedure for eacn run was to vary the membrane airflow so 
that one oktained both positive and negative pressure 
gradients at the chamber orifice. This provided a =t«:ore 
precise location of tne proper entrained flow rate. It is 
of interest to note that, due to hardware limitations and 
mre targe ancunt of air entrained by the jet at the Lurthes 
axial Ilccations measured, nozzle exit velocities rere 
dimited to around 100 feet per second at the further axial 
mecations. Tatss also ~aterected the  hunber OL useful 
locations for whistler measurements, aS a strong whistler 
oscillaticn required a relatively large exit velocity, say 


300 feet per second. 


Waa heduction followed standard procedures for the 
Memturt and oritice. These are presented in appendix C. 
Entrainment was plotted as M/MO(total mass flow divided py 
the primary mass flow), against X/D, axial location divided 
PyenozzZie diameter. forking curves of chamber orifice 
Peceemrce gradient are presented as pressure daftference 


against Hy“0, in the example plotted in figure 10. 
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Ce. FLOW VISUALIZATION 


weg N2oualization of the nozzle jet lip interaction was 
performed in two ways: high-speed motion photography of a 
meegeOr tufts at the nozzle lip, for the purpose of 
discerning a pattern to the interactions; and a mixture of 
lamp black end kerosene inside the nozzle collar in order to 
observe the behavior of the separation reattachment 


Mmocation. 


Photcgraphic equipment consisted of a Honeywell Fentax 
Single reflex camera, a Bolex Supreme handheld 16mm _ rnovie 
Camera, and a xedlake Hycam 16mm high.~speed novie camera. 
ietesclte Ling WaS Simpiy constructed of cotton threads held 
by a masking tape ring, and a photo of the ring is presented 
in figure it. Observation in the lampblack experiment was 
performed with the naked eye and no peripheral equipment was 


PececscaLly. 


Photograrhic tests were to pe done in three phases. 
High-speed still shots were first taken to determine the 
usefulness of the tufts for visualization, and to obtain the 
approximate reguired tilm speed for motion photography. 
Next, medium speed motion pictures, 64 frames per second, 
were taken to confirm the necessity of the high-— speed 
camera. Firally, high-speed motion pictures were taken to 
weotaily ccntrol the tuft motion such that userul 
Oobservaticns could be made. No particular collar settings 
Or exit velocities were used in photographic work. The only 
reguirement was to produce the oscillation phenomenon, which 
Peemmever GQifiLicult. AS it turned out, a lack of tine 
prevented the completion of the high-speed photography 


phase. 


Sy 





Lanpblack tests were performed in two ways. The inside 
of the ccilar extension was liberally coated with the 
lampblack mixture, and either exit velocity varied at 
@onsctanht CCM cetting, Or coliapwisetting varied at 
constant exit velocity. Observations were all visual, none 
quantitative, save noting the ranye of exit velocities 
covered, which was from about 100 to 300 feet per second. 
The lampblack mixture was observed to pool as illustrated in 
figure 12, and the downstream—most edge or the pool was 
Miter preted damtmda2cating the Jocatton on tree collar of 


reattachment and Separation. 
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Nozzie thrust efficiency results ate presented in figure 
13, plotted for various collar positions. The whistler 


oscillaticn was produced at each collar Setting for “arc 


» 


Fan 


whole range of Reynolds numbers’ shown, although sudden 
meequency JUMPS did occur in this range. These freguency 
juaps did net appear to change the nozzle thrust or thrust 
efficiency in any discontinuous manner. Instead the trena 
WaS as indicated, progressively decreasing nozzle efficiency 
With increasing collar extension. Whe zero collar extension 
efficiencies were found to oe about 6 percent lower than 


those found in a previous report (reference 8). This ais 
believed due to the extremely simple, unflared configuration 
Seco mOzZic, which may be seen in figure 4. reference 8 
nozzles were divergent nozzles and as such should be more 
ereicient. Important to note is iva the whistler 
OSCliilations tay ve obtained without wuch reduction in the 


Mascrc nozzle efficiency. 


Bntrainment results are presented in figure 14, plotted 


together with the results of Hiil and Greene (reference 1) 


mon COM CaLiscn. AS may be seen, in the case of the basic 
jet agreement is sood guakitatively pure poor 
Simemecatatively. In the case of the whistier nozzle 


agreement is not good in any way. It is believed that the 
Foor guantitative comparison vas due to interference of the 
deep orifice cylinders with the jet. These create a 
pressure gradient from the orifice to the chamber exit which 
guite possibly affects the readings. As noted in the 
experiment Seoction, “and Sin *hppendax’ C,. these “orifice 


cylinders constitute a major difference in configuration 


She 





from the chambers of reference 7 and 8. Additionally, in 
the case ci the whistier nozzle, it is believed the method 
of enciocsing the jet up to the axial position of measurement 
May interfere with the sound - jet interaction (especially 
the cloth porous cylinder) andthe formation of a larde-~ 
Pease tCurdulent structure. Figive 5 i Jerr ates setae 
readings taken were relatively constant with Reynoids 
Meer. AS Noted by Ricou and Spalding, the entraiument 
ratio can be a strong function of Reynolds number in certain 
ranges of Reynolds number. It is interesting to note that 
this behavicr 1S supported by the present notion of jet 
turbulence structure and its development as exit velocity 
increases, hich postulates that the matior contributcr to 
Pm@eraoanment, the large scale structure, dissappearts as 


Reynolds numter increases. 
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Lack of time precluded testing at more jet axial 
locations, as each location required a different orifice 
Gylinder (see appendix A). Of the originai folic Orifice 
cylinders constructed using reference 8 data, only two were 
useadle as the Spreading angles encountered averaged about 
28 degrees, cr 10 degrees less than those used in reference 
ore Himeroraders to be effective, the orifices had to re 
positioned at axial distances further than they were 


designed fcr. One of those was modified to obtain readings 


miele s GSSign position, also. The larger two were never 
erfective. Ahnotner problem encountered was the mass flox 
Memstaticn cr the porous cylinder air Spo This 


mestetcted the useable nozzle flow at axial locations kEeyone 
mee ot  ahout 6 to such low values that a whistler 
Seciitation strong enough to affect the entrainment could 


not be vreduced. 


He was hoped that high speed motion pictures t0ould be 
taken of the tuft experiment; however, this phase of the 
photographic experiments was not reached. The effectiveness 
Moemeere tufts a2n iiandicating the circumferential colla 
positions of attached and detached flow was proven in tne 
high-speed still picture paase. Motion pictures with a oi 
frame per second camera showed that a much higher speed 
camera waS necessary to Slow the motion such that  frotion 
pictures would be useful in slowing the tuft motion to a 
reasonable speed. The necessary film speed was estimated as 
about 800 frames per second; however, the pictures were never 
taken. Scme representative still pictures are presented in 


Ergure 16. 


Lampblack tests appeared to be productive aS indicators 
of the general behavior of the collar attachment point. 
fdicattions are that a sort of Coanda effect is involved, 


memmered by the fact that while the jet 1S attached at any 


4 





point on the collar there are other points where detached 
meow altowS a backflow into the collar jet gap,which tends 
to relieve the Coanda effect causing a cixcumferential 
Shift in the point (or points) of attachment. At a fixed 
collaz extension, with increasing jet velocity, the attachment 
point was ohserved to move upstream to a limiting position, 
just as the Coanda jet described in reference 10. However, 
Temi kea CXiMMcHocity and inckeasamgecoijlar extension a 
different behavior was observed. The attachment point 
seemed to move With the collar until a point was reacned at 
Which the frequency of oscillation jumped. At that instant 
the point of attachment also appeared to jump suddenly hack 


to its original position relative to the jet exit. This 


ay 


ehavior indicated some interdependence of the whistler 


attachment phenomenon and the oscillation phenonenon. 
In addition to the above experiments, two other 
Observations were nade which verified those of Hill and 


Greene. 


Parst, it ue emeoted= tiharemextencioemmer the whistier 


collar seemed to increase mass flow througn the nozzle. 
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This was confirmed with static pressure measurements in the 
mozzZie "throat'. These indicated that tne collar was acting 
as an overexpanded diffuser, resulting in an effective exit 
area Slightly greater than the ‘throat' cross secticn but 
Meas than the collar cross section, and further resulting in 


a greater mass flow througn the nozzle. 


Second, it was noted that a correspondence betwee2n 
mecenuency Of oscillation, jet exit velocity, and collar 
extension seemed to exist. It was found that by gradually 
extending the collar as plenum pressure wWaS increased, one 
Souter *tOliow’ a frequency or oscillation. RECO EC imine. 
collar extenSions were measured and jet velocity caiculated 
at fixed frequencies of oscillation, and it was found that 
these fixed freyuencies corresponded to constant ratios of 
collar extension to jet velocity, thus confirming the organ 
Pee character of the whistler oscillation noted by Hill and 


Greene. 
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vit CONCTUSZONS  Aaip RECOM SRATIONS 
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The thrust efficiency and mass entraingent. 
characteristics qf the whistler nozzle have been measured, 
in the case cf the entrainment, using a technique neretofore 
moeesesei On ah unsteady jet. Further, observations have 
weemenade Which shculd contribute to a further understanding 
of the whistler flcw mechanisn. 


Severdi directions of further work are indicated 
First, it is recommended that a thrust measurement device 
uSing @ pendulum arrangement, Similar to that used in 
reference to he constructed £ OG a nore accurate 
determination of the whistier efficiency. A device of this 


type might avoid the hysteresis problems encountered with 


Mower tieUusSstS cn a conventional thrust bed. Lt “1a Ssaeaisc 
reconmended taac dad “WhasSstler omozzlie of Smaller” finroat 
diemeter be constructed 2e permit nass entrainnent 


measurements , at further jet axial locations, in the 
entrainment chamber, at the high exit velocities required to 
Peomdwce a Strong oscillation. This will also permit further 
study of the possibiiity oi whistier measurements in the 
chamber cver a greater range of Reynolds number. Next it is 
recommended that further flow visualization experiments be 
conducted. Specifically, high speed motion pictures of the 
tufts used here, and high speed scniieren photography cf the 
jeterswould prove wseful. Finally, it is reconmended that 
other unsteady nozzles be tested in the chamber, as the poor 
results oktained here with the whistler are only preliminary 
and may not ke indicative of the usefulness of the chanber 


Mrenenon-acoustic efrects such as swirl. 
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VIIT. INTRODUCTION- THE RING WING 


does prescnt work concerns a hollow cylinder,  ameially 
aligned, in supersonic flow and undergoing small amplitude 
Secitiations in angle of attack. The intent is to calculate 
the resulting pressure distributions and aerodynamic mcnen 


on what might be mcre properly termed a ring wing. 


Previous inveStigations, i.e. references 12 and 13, have 
eypizes the method of characteristics towards the separate 
Seweculecticn of the inner and outer flowfields, as tell as 
pressure distributions and generalized aerodynamic fcrces, 
#Or such a cylinder in supersonic flow whose walls are 
undergoing small amplitude Sinusoidal oscilijations, or panel 
flutter. These previous investigations resulted in computer 


programs to cerform the necessary caiculations. 


The primary thrust of the present work was the 
modification of the previously created programs and then the 
Sewpinadticn Of theiz resuits to calctuiate the unsteady 
pressure distribution, lift and wmoment on the cylincer 
undergoing angle cf attack osciilations. Finally, using tne 
moditied programs, a parametric study of these items was 


conducted. 
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IX. PROBLEM FORMULATION 


a oo ee ow SE owe 


This development is ciosely based on several previous 


papers , 1.€. references 12 and 13. 


Soneidem a carcular cylindrical Shell, ih -supersenic, 


inviscid, adiabatic flow, whose axis is aligned witha the 


freestream, and whicl is undergoing small anplituaude 
Weer iations in angie of attack, about the Z axis. {figure 
a) 


MeemeGover ning equation for this £flogm Situation 1s tue 
Jinearized unsteady potential equation, in ey bindiiea: 


Geordinates. 


Oieecerst O., t+ 1 Oo +1 G5 - 2M OL. 7 1 oo = 8 


fee relevant bcundeaery conditions are the flow cengency 
Sond1tion and Sommerfeid's radiation condition, 1.€.,thnat 
ee 


disturbances will propaga away from their source. 
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One further condition is necessary,along the X-axis, which 


amounts to setting the pressure disturbkance at the axis 


egual to zero. 


The linearized form of the flow tangency condition is: 


Which must be satisfied on the cylinder inner and outer 


Sue tace. 


AS Statec, the boundary condition at the axis iS a zero 


pressure disturbance: 


The linearized form of the pressure coefficient is as 


fer jlows: 


uSing cylinder iength and 


Arter nondinensionalization, 
the above equations 


length divided by freestream velocity, 


become: 
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(1-M?) + + 1 + 1 - 2m? - M? 
ee ee ge ot Oe ee = 
i ie 
= 2 + 
: a ee 
Sub ject to: 
¢-| = ah+ ah js xl 
r=R ) ax 
One now assumes a cylinder oscillation of the form: 


Peewee) — Z2{(x) cos (6) cee 


where the nondimensionai fredquency, kK, is: 


Mimemeet lection amplitude, Z2({x), for the case wuniher 


Senpsigeration, 31S of the form: 


mapa = et x = x) 
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Following the assumed form of the deflection, the 


perturbation potential must be of the forn: 


fee, 0,t) = 6(x,r) cos (@) Ee 


and, after substitution, the basic equations hbecone: 


—~M? + + 1 + (k*M*-1 - i2kM =O 
roi) oy Ca = Y (} : ) 9 On 
Coem 7 Ox tek? ) 

emp yect to: 
>| eee pO eax SL 
a 
r=R 
and: 
ea) 0 
P y=0 
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X. METHOD OF CHARACTERISTICS 


LOL SuUrerSoOnic Flow, the linearized perturkation 
Meysent2zal equation is hyperbolic, and has characteristics 


Pimen Satisfy the following differential eguation: 


(Me-1) d*r - d*x = 0 


Since tte cnaracteristics in supersonic flow correspond 
to the tiach lines, then for ds, the differential arc length 


mepong a characteristic: 


ds M as M 
An ar Sitrary HENCE LON ee © (x, fo ew rias , these 2c browong 


derivatives alcng the characteristics: 


a = Ee ax + i eS /M*—-1 ee dl ~ 
as A ds M M 


Defining dst and ds2 to be the differentiai are lengths 


mmogetne tert,and rigat-runking Mach lines respectavely, 


then: 
ae Ol ' Hy cle 
-dsy ds, 


a2 





The cross derivatives become: 


SOlvingd mow for Fxeand Fre 


1 ee M (F_ + FF) 
oye 86 é 


Now the basic equations may be written in terms of b 


Mmamwies Gerivatives along the Cnaavacteristics as follows: 





= =v - Te eer el - ikM + o_) 
eg a ee jf 7 ae, 
(elo. =— >.) = 24. + 1kZ 
r 5 1 2 
Cl = -2(¢, + ikd) 
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biese cQuations have been put into finite difference 
HOrm and frogrammed in Fortran IV for the separate cals@s o 
exterior flow and interior flow of a cylinder whose walls 
Were undergoing small amplitude sinusoidal oscillations. I 
the interior flow program the X axis boundary conditicn i 
Perse iod Gi a sittall Gnner cylinder in place of the ak#s, to 
smrcumvent problen= encountered should r go to zero. FOr 
detailed acccunts of the equations and programs Rotee 20 


EFeferences 12 and 13. 


The primary worx here was the transforsmaticn cL the 
aforementioned programs to the case under consideration , 


H5€., a cylinder undergoing small amplitude oscillations in 
plgte Of attack, and then a parametric variation study of 
the pressure distributions and generalized a 


force. The programs were nodified such that all 


rg 


= 
Cc 

Mete swith an angie of attack amplitude of 0.1 Yadiens. 
Bice the baSic egquation iS in a linearised form, the 
conversion of any  ©esuits to “Other. angle. On fecrac" 
meet udes iS a Simple matter “or multiplication. For those 
interested, the programs resulting from chahyes reuyuired to 


the reference 12 and 13 programs are iisted in appendix D. 
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AS Stated above, the 


coefficient amplitude ji 


aerodynamic force, One. 
follows: 
meee | OC Vents ax 
e os Pm ~ 
wee = © Ax WV. (x) 
a 


Where Crm is the pressure 
from 
deflection node, 


Geer tection mode. 


the m'th axial deflection mode, 


anc A is the amplitude of 


Posi Vesecoeud veda we ae 


ASLL UD IeLON Sando eh 


the 


pressure 


generalized 


Tn its most general forn, QOmr 1s 4s 


coefficient amplitude 
Psi 1S tiie 


the 


O@eta~es the following forn: 


Cry tx) dx 
fe) poe 


oa 
2 


4th 


Fe eG 
2 ie Seed 


axial 


For the special case uncer consideration, 





Thus , in the case under study Q has the special 
interpretaticn of being the moment coefficient amplitude 
about the peint Xc multiplied by the angle of attack 


mires tude, daipha, which is equal to 0.1 radians in all cases 


Peesented . 





several cases cf comparison were considered towards the 
evaluation cif program results. Figure 18 presents two cases 
of compariscn for pressure distributions. The outer radius- 
meee encth ratio was taken to such a value (10) as to 
approximate a Be el yh plate se Korey LAE The Luner 
Beadiwc-~ctoO-lencth ratio was taken to such a value (9.75) AS 
Memecalse da reflection on the inner surrkace of the outer 


eyilander. 


meee suit icivently large “radius =to-length ratzos tie 
exterior pressure distribution should approxinate the 
[PmcOhetical value for a flat plate in supersonic flow, i.e. 
2o/f8- Figuce 18A presents this comparison, and as can be 


seen agreement is good. 


For sufficiently iarge radius-—-to-length ratics the 
moter ior pressure GIStEL Outen  siould yea pprOci nd se eee 
theoretical value for a flat plate. de OT a) 1) ae roe 
Mier cylinder iS Gk such a radius-to-lenyth ratio as to 
GCauce a Mach wave to reflect upon the outer cylinder inner 
surface, a flow situation exists which should approximate a 
flat plate ina free jet, where the mach waves are reflected 
maotmeeecne jet pyoundary. This is due to the. boundary 
eondition imposed on the inner cylinder, Cp = Q. Reference 
Home has treated this case in analytic form for smail 
frequencies. Figure 18, A and B, presents a comparison of 
Meegream results and the results of fererence 16. AS Can be 
seen agreement is good. | 
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Reference 14 has treated the case of a nonascillating 
ring Wing in supersonic flow. Figure 19 presents a 
comparison of prodqram results and those of reference i4, for 


the case shown. As can be seen agreement is fair. 


In veference 12, which develoved the hasic inner fiow 
field program used in this report, the occurrence of a 
numerical instability in the pressure distributions, under 
certain conditions, was noted. Tt was found ethat seis 
@isturbance, which occurred only after a reflection, was cue 
to an undefined constant in the subroutine which compntec 
the boundary condition at the axis. After defining the 
Memecams. Correctly time pressure distributions were Seen to 
settle down quite nicely. Figure 20 presents a before ana 
after conparison of pressure diStrrbur venus OS a 


representative case. 
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Figures 21 through 23 present various parametric studies 
Sistine MOnCH< Ccoeificient, 0, fer cases both with and 
Without a reflection on the inner cylinder surface. In 
these figures the coefficient plotted 1s actually the 
a@ifference ketween that calculated for the inner cylinder 
Surface and that calculated for the outer cylinder surface, 
divided by seie- “angle sof S.attack sanplutidea. Thus the 
coefficients plotted are Ota! monent and pressure 
coé€ificients, as opposed to those plotted earlier, which are 
the coefficients for the inner and outer cases separately. 
AS can be seen from figure 22,for a center of rotation at 
the leading edge the cylinder has definite geometric 
regions cit stability with respect to angle of attack. 
pure 21 1llustrates the effect of increased frequency on 
the moment coefficient. EUGUGe 23a iet rates aero tea 
pressure Gistribution including one reflection, from which 
emewean infer the reason for the behavior of Q in fiqgute 22, 
and also the. behavior of Q with a varying center of 


Botatlon. 
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XIII, CONCLUSIONS AND RECOMMENDATTONS 


— weet OEP ar eee ee ee See ee ow ae = gene TES TUR Cee TUS Tee ewe qeeme See Oe OE Oe ee 


The computer programs of references 72 and 13 have heen 
PuiceessS*ully adapted to the case of a ring wind me 
Supersonic flow undergoing oscillations in anale of attac’. 
Available cases for comparison have shown good agreement, 
Meecating the validity of the method. The programs rake 


possible the relatively easy calculation of the unsteatv 


[eres and moments of a rind Wing in supersonic flow, tor 


Sten aDdDlications as tubular artillery shells (S127 
spinning e@oular projective. ali nougiee spilge tay smc. 
aerodynamic effects not accounted for in the metho. 


Further evaluation and parametric studies are reconmenie: 


fhoeenea More complete verification of the computer programs. 


66 





APPENDIX A 


ASAT ROPY Oo THRUS’ CALCULA Tags 


Me-2-TOpie thrust calculations made use of Sitanrmara 
isentrovic relations. Calculations vere hase@ on  pdlepum 
Peeo>tre and nozzle exit area. The formula for isentropic 
imomtSst 1S arrived at as folldes, Wocele thrust Way eve 
Mertten in the following forn: 


T = wiV = (o0AM°YRT 
S ex SS —~ 


where Wis nozzle mass flow, Ais exit area, and one has 


har bd 


X 
@eetred the speed of sound relation. Exit Mach nuzber may 


be written in terms of pressure ratio as follows: 


ic 
I 


eee ae) 


ex yTh ae 


Applving this and the state equation to the formula for 


ms One obtarans: 


agit 
ID eee 2 ee =a) 
Now, once one has measured plenum and ambient pressure, 
and thrust, Texp, one can calculate the isentropic thrust 


efficiency: 
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Re PENI xX 5 


PE Rete Oka Po eh LeU ere 5 


Nass flows and entrainment ratios were calculated 
accordinzjy to the equations and procedures sat out in 
reference 11, the ASME Power Test Contes. The equation for 


mass floyv, from reference 11 is: 


w = 359 c f d* Fa Y Vvho ikoy Abbe 


rh 


Nye B=d/D 
dé and D d | D 
ue ee 


"a Thermal Expansion Factor 
Wiwausor a Venturi) = Adiabatic wiaxpansi on actor 


feeckine pressure Jirop acioss the @evice, in inehes oT 
water 
9 Density of the fluid being metered, pounds per 


Euigic foot 


TF wt is the entrained flow and w2 the primary, then (#1 


fheesua Venturi, W2 an orifice): 


=I = 


= A vh 
meee 7 (> See 
; cea) as Fa5 yt hp 
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on-OComethic factors For the primary flow; 
dischargqe coefficient, 


such AS 


could be reasonably approximated for 
all tests as: 


C2=.65 


| 
py 
th 
I 
a 
e 
>) 


Y2=.99 


In the case of the entrained flow these factors had to 
be changed for each test axial location, due to the widely 
varying flow requirements. The aqeometric factors above had 
the followino fired values: 


Mite JVo28 dai=.75 P= Ve0yoe G21. 1% 


Yow, entrainment ratio is commonly written as: 


fia wv. + WwW. = 


W 
a So 5 


0 2 


tes 


teang the previous formula for “wl/w2, the ~entrainnen 


ck 


ratio may now be written as: 
M =AvVh, +1 
ae W 
Meee A as the combination of the ratics of fac Toms 
previously listed. This equation vas used to calculate the 


foze oo eontrasanment ratios. 





hee egy tC 


ENTRATNEINT CHAMBER DESIGN 


The entrainment chanber design was closely hased on 
Peer OuS Fork by Rico an@ Spalding (reference 7), “ana 
Peschke (referenca 8). Peschke'ts report especially was a 
qood sonrce of specific design data, while Ricou an 


Sealdingd was primarily a source of the desiaqn Principles. 


The entrainment chamber consisted essentially of a 


na 


mS Wwtera! woreus cylindar, closed at one end with a solid wal 


through which the nozzle protruded, the other end being o 


eS) 


2, 
=) 


to the outside. Orifice plates of svecific dianeters ve 


ry! 


° 
constructed and mated to cylinders of specific length, 
enge h and jVeteter ) Speci n red inccor. Ings) tom Cie madara 
available in Peschke. These orifice cvlinders were of such 
outside diameter that they fit closelv inside the vorous 
cylinder, into which thay were inserted for the measurement 
of mass flow at the Veo. Fele, axia le posi frons 
corresvonding to the orifice diameters. The purpose of the 
orifice was to seal the open end such that entrained alr 
could only enter through the porous cylinder. The purpose 
of the orifice cylinder was to seal off that nortion of the 
porous cylinder beyond the jet axial position where flow was 
being measured, The porous cylinder was surrounded 
everywhere, except at the open end, by a box, also encasina 
the nozzle plenun, which served as the nembrane air supply. 
The entrainment chambsr was about 4 by 4 feet square and 


Merce high, Stipported on 3 foot leqs. The porous cylinier 
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was 30 inches in diameter and 25 inches long. Orifice 
diameters used, corresponding to the Listed fet axial 
positions, were as follows: 
Na aneter x7) 
eas IS: 
2.0 oa) 
Oa) i268 
Paemsentrainnent chamber box; legs, orifice Glates amd 
Me-ous cylinder frame work were all constructed of weod. 
mee) OLitTice Cylinders wer2 made of fiberglass sheet and the 
Merous covering to the porous cylinder frame “work was two 


meyers of dacron crepe cloth, as suqgested in Peschke's 


report. 


Q 


— 


One Critical chamber design parameter vas orifi 
Giameter. It is clear that too small an orifice vould 
interfere with the jet being measured, while too larga an 
SureLcesfouldgnot proparly restrict the entrained Flow to 
entering only through the membrane, thus preventing the 
creation of a readable pressure gradiant across the orifice. 
The approach to this proptem in both references 7 and 6 was 
to find an optimal orifice size at a specific jet axial 
location hy testing several orifice sizes. fo minimize 
Orifice construction, this investication took a diftferent 
approach. That was to use the reference 8 orifice data as a 
Sesto moe pOlnte anidmvaly thiemori fice “axial location, ~ such 
that a readable pressure gradient was created while not 
Interfering with the jet. In both approaches the attempt 1s 


to optimize the orifice diameter, jet axial location 


combination. 
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One primary difference between this chamber's geometry 
anc that of others (references 7 and 8) should he noted. As 
Shown in figure 8, exit orifices were positioned in the 
porous cylinder at various distances from the nozzle. That 
portion of the porous cylinder extending above the oritice 
plate was sealed by a fiberglass cylinder attaches} to the 
@eaftice. PIMs) thas Sivuet lon, tliat port on meet neta. 
beyond the orifice, and inside the porous cylinder exit, was 
mot in a truly free jet flow situation. This was as opnosed 
aomere reference 7 and ® configurations vhere the “exit 
Meperce Corresponded to the exit of the porois cylinder. 


Thus in those experiments the jet past the orific= was in a 


meudy  rree jet flow situation. The configuration used in 
mois renort was adopted in order to Facilitate the 


meagouberent of mass flows at many ayjal locations. It sas 
mepea that this chamber's qgeonetry would net affect The 
Chanber readings oareatly, although sone effect was to be 


— 


expecten, 


An inportant probtlen, again vith regard to the chamber 
orifices, concerned the measuremen of the pressure 
differential across the orifice. Roth reference 7 and &@ 


indicated that a measuremant accuracy of .0901 inches of 


Water waS required. Reference 7 used a micro manomete 


H 


while reference 8 used a pressure transducer for thi 


P) 


purpose, After the entrainment chamber was constructed 
experience shoved that stitch accuracies were required. aes 
investigation used an alchohol micro manometer which 
provided an accuracy of .0901 inches of manonetric fluid. 
This accuracy was obtained through the conbination of an 
inclined tube, hairline indicator, magnifying lens and a 


micrometer adjustment to indicator level. 
Flow measurement devices for the primary and entrained 


Ellow, were an orifice and venturi respectively. These were 


mayor eucted and used according, to WASN= Power I[=sst €oadges 


eZ 





specifications (reference ee Rased on these 
specifications, an inclineli water manometer was found to he 
necessarv for measuring the orifice pressure drop, while a 


vertical mercury manemeter was sufficient for the venturi. 
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bDP TT DED 


PROGRAMS LISTINGS AND DOCUMENTATION 


This appendix contains basic Progra (“ard Si nous 
Pescriotions for the two programs used, from references 12 
and 13. The programs also are listet in their nod@dified form 
for the computation of a ring wing at angle of attack. The 


proqran statenents relevant CO> ths) modi teas ong 


tae oo 


indicated in the program listings. Also indicated in the 
interior flowfield program is the statement whose omission 


caused the numerical instability noted in reference 12. 


The inner flow field program consists of & Subroutines, 
and a main program which manipulates each of these. The 
action of each particular subroutine is indicated in the 
program listing. Inputs are made by way of a Fortran 


Namelist. The input parameters are, in prover orier: 
FSTRNN Free stream mach number. 


mee nO The reduced frequency. 


RO [iVcmouTenvcvlinder TadxusssOutength rat io. 
! 
Ral Hie Atiner CVliInger tbadius to Length ratio- 
NGRDFN Grid fineness, or the number of gtid points 


Paleeneon the first Light running mach line. 


N Circumferential node number. 
4 Axial mode number. 
M7 The m in Omn, the generalized aerodynamic force. 
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TPRINT An output paraneter. 0 indicates only basic 
output. 1 indicates output including such things as Phi and 


its derivatives at each flow field point. 


Of the above inputs, in the modified program, “he tnputs 
N, “4, and M1 must be set to 1, as the case under 
consideration is not panel flutter but a cylinder at anale 


of attack. fhe other inputs are arpitrary. 


15 





Ee EEEEECECECEEECEEEEECC ECL CECE LECCE CEC EL ECCCE CCC CCCEEE CECE 
iS CYLINORICAL SHELL/PANEL FLUTTER PROGRAM ; 
CEC CECCCCEC EE CEECEECECE CC CC CCCCCCCCCCCCCCCEC CE CECE CEC CECE 
THIS IS THE MAIN PROGRAM 
COMPLEX PHICS00) sDPHIS1(409) sDPHIS 21400} 
DIMENSION X{400} »R(490},XPT(400) sCPRI(400),CP1(400) 
DIMENSTGN DUMMY (400) ,DUNCE (400) 
CCMMON/ BECKS FSTRNN, DELTAS » DX, DH» OSTSTR RO ,RIyREDERO 
COMMON/BLCK2/ Al, A2,A2¢h4sRNyRMeRML FACT Ly FACT 2 
COMMON/BLCK2/ NGROFNSNPTS,LCOCUNT yg IHAVEP »NyMyhMl,I PRINT 
CCVMON/BLCK4/ PHIiyDPHIS1] »OPHIS2 
COVMMON/BLCKS/ XyReXPTeCPRICPI 
; UATA PI/3.141593/ 
_ ENON SU ee eC oar Van) ee Oe Ao 7 ree 
C 
Bio = beck 04 
LER = 0 
CALL INPUTI(T ER} 
en el EE eee 2 GO LTO. 2 
aS CONTINUE 
IHAVEP = j 
C POEUN TF =" 4 
: CALL INITAL 
C 1 CONTINUE 
Gee COMPAR 
C Poin on lO POINT GNOTHE INITIAL RIGHT RUNNING MACH 
: LINE? 
Peete ECGUNT .EG. 1 ) .ANDe. CIRAVEP «LE« NGRDFN )) 
3 Ht GU 1.0. 
: Womiitjeeor lb POINT ON THe CUTER CYLINDER SURFACE? 
. MaotelitaAver sE0%° 0 } GQ 10 %& 
A Pies o nt OtNiON THE INNER CYLINDER SURFACE? 
A PoetetaAVEP Set O. NGRDENG? GO 1G 6 
5 Poe THe GRID POINT IlSetN VHeSGeNersAt FLOW FIELD. 
ewer GENFPROILER) 
Pe lek Pee Ose lL. p@Go. 10 12 
TAVEP = THAVEP + I 
2 Gear | 
Beal MACHLN 
Petey Ee? stO. NGRDEN ) "GO t0e4 
THAVEP = IHAVEP +t 1 
2 G@ 70 1 
C THIS IS THe LAST POINT ON THE INITIAL MACH LINE. 
c 
4 tCCUNT = LCOCUNT + Il 
Paver = 0 
, Gist 4) 
: Sy GS eA 
C oath S GRID POlNTe PAST Ve eEex)eLE CYLINDER LENGTH? 
£ 
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WY 
N 


FSTRMN;s REDFRO,RO,RISN 
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WRITE(6;114) 
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eo CIRPOUM-ERENTIAL MODE *, 
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CCMPLEX DPHIDR,DPHIDX 
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The outer flowfield program, from reference 13, consists 


of 5 subroutines and a main program to manipulate the 


1} 


subroutines. The action of each particular subroutine 
indicated in the program listing. Inputs required to the 


program must be on two cards, and are as follows: 
DATE 
MACH Mech number. 


RADTUS Padius to length ratio. 


ie Reduced freguency. 

MFREO Axial mode number, 

NR The r in QOmr, the generalized aerodynamic force. 
NFREO Circumferential mode number. 

FINGRD Grid fireness, essentially the number of arid 


mes a vonGg tre cylinder surface. 
$ y 


The run date is typed onto the first card. The second card 
contains the above ovarameters in the followine format: 
Pee pet. 3, 415. OF the above parameters, in =the 
Were 12¢ “program, MFRFO, “NR, end NFREQ must he set to 1, as 
the case under consideration is not panel flutter but a 

ther 


cylinder at anagle Oe a tia 0 paraneters are 
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